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Abstract 

To clarify the yielding mechanism of small 
hydrocarbon molecules in chemical sput- 
tering between hydrogen and graphene 
sheets, -we made classical molecular dy- 
namics simulation with modified Bren- 
ner's REBO potential which we proposed 
to deal with chemical reaction. As the 
simulation model, we prepared more real- 
istic physical system, which is composed 
of 160 incident hydrogen atoms and ten 
graphene multilayers, than our previous 
model. From the present work, we found 
the following fact: breaking the covalent 
bonds between carbon atoms by hydrogen 
does not play an important role during de- 
struction process of graphene structure, 
but momentum transfer from incident 
hydrogen to graphene causes to destroy 
graphene structure. Moreover, it is found 
that almost all fragments of graphene 
sheets form chain-shaped molecules, and 
that yielded hydrocarbon molecules are 
composed of carbon chain and single 
hydrogen-atom. 
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1 Introduction 

Plasma-carbon interaction yields small hydro- 
carbon molecules on divertor region of a nuclear 
fusion device [1] . Diffusing from divertor region 
to core plasma region of fusion device, gener- 
ated hydrocarbon takes energy from the core 
plasma. Therefore, the hydrocarbon is some- 
times called "dust" by the nuclear plasma re- 
searchers. To remove an ill effect of the hydro- 
carbon on core plasma, ■we need to kno^w a cre- 
ation mechanism of the hydrocarbons. Ho-wever, 
the creation mechanism of the hydrocarbons has 
not been elucidated yet. 

As the first step to clarify the creation mech- 
anism, we investigated, by computer simula- 
tion, collision process of hydrogen atoms and 
one graphene sheet, ■which is regarded as one 



of basic processes of complex plasma-carbon in- 
teraction in the previous ■work [2]. In the pre- 
vious simulation, we used 'classical' molecular 
dynamics (CMD) algorithm ■with modified Bren- 
ner's REBO potential which we proposed to deal 
with chemical reaction between hydrogen and 
graphene [2,3]. From the previous work [2] in 
which an incident hydrogen kinetic energy Ei is 
less than 100 eV to compare with experiments 
of divertor, it was found that an hydrogen- 
absorption rate of one graphene sheet depends 
on the incident hydrogen energy, and that the 
collision mechanism between a graphene and a 
hydrogen can be classified into three types of 
processes: absorption process for Ei < 5 eV, 
reflection process for 5 eV < £"1 < 50 eV, and 
penetration process for Ei > 50 eV. Moreover, 
it was also found that when hydrogen atom is 
absorbed by graphene, the nearest carbon atom 
overhangs from the graphene which we called 
"overhang structure" . 

Based on the above results, simulation model 
are extended from a single graphene sheet to 
multilayer graphene sheets in the present work. 
We adopt three cases (5 eV, 15 eV and 100 
eV) as the incident energies of the hydrogen 
atoms El. Each value corresponds to the ab- 
sorption, the penetration, and the penetration 
process of the collision between single graphene 
sheet and hydrogen atoms, respectively. More- 
over, we use the same simulation algorithm and 
the same interaction potential as the previous 
single-graphene sheet simulation. The aim of the 
extension of simulation model is that we reveal 
a more realistic sputtering process of graphene 
sheets and hydrogen atoms than the previous 
work. 

In §. 121 we review the simulation method in- 
cluding the interaction potential which was pro- 
posed in the previous work. We discuss the sim- 
ulation results in §. |21and summarize a conclu- 
sion in the last section §. ^ 

2 Simulation Method 

The simulation algorithm and the interaction 
potential are the almost same as the previous 
work [2]. First, we review the interaction poten- 



tial among hydrogen and graphite atoms. The 
original Brenner's REBO potential [3], which 
was proposed on the basis of past simulations 
[4-6] , has the following form: 

(2.1) 

where is the distance between atoms i and 
j, is an attractive term, is a repulsive 
term and the function bij{{r}) includes all ef- 
fects of molecular orbitals. However, if chemical 
reaction occurs, the REBO potential breaks en- 
ergy conservation. Therefore, we proposed the 
following new functions expressing conjugation 
effects [2]: 

carbon 
carbo7i 

+ E rir,i)CN{Nl), (2.2) 

where f is a cut-off function for a distance be- 
tween atoms, and 

Cn{x) 

( 1 ifx<2, 

i+cos{^(x-2)} if2< a; < 3,(2.3) 

[ ifx>3, 

nl = E rK)-rKO 

= E /'(^fc^)- (2-4) 

Modified Brenner's REBO potential conserves 
the total energy during the chemical reaction. 

Next, we explain the simulation model in 
the present work. We prepare ten graphene 
layers each of which is composed by 160 car- 
bon atoms. All the layers are perpendicular 
to z-axis and have periodic boundary condi- 
tions for X and y directions. Distance between 
graphene sheets was set to 3. 348 A in the initial 
state. Graphite atoms at four corners of the first 
graphene sheet from the top are fixed in all simu- 
lations. In the second graphene sheet, the center 
graphene atom is fixed. After the third graphene 
sheet and later, the above fixed condition is re- 
peated in all layers. As the initial condition of 
the graphene layers, we prepare the equilibrium 
state of carbon atoms in all graphcnes with a 
temperature 300K. 

We shoot one incident hydrogen atom at a 
time every 5 x 10~^^ seconds from z — 120A 



plane to the first graphene which is located at 
z = 3.348Ax(4+ 1/2) = 15.066A, until the to- 
tal number of the hydrogen atoms becomes 160. 
We make the three simulations each incident hy- 
drogen of which is set to 5 eV, 15 eV or 100 eV, 
respectively. The x-y coordinates of the incident 
hydrogen atom at z = 120 A plane are given as 
a pair of random numbers ever time. Each in- 
jection angle of the incident hydrogen atom be- 
tween z-axis is set to a uniform random number 
from 0° to 60°. 

If an incident hydrogen atom is reflected by 
the graphenes and reaches z — 120A plane, the 
reflected hydrogen atom is omitted from the sim- 
ulation. If the incident atom penetrates the 
graphene sheets and reaches z — — 120A plane, 
the penetrated hydrogen atom is also omitted 
from the simulation. 

To integrate the equation of the motion, the 
second-order symplectic integration method is 
used in our simulation [7-10]. 

3 Simulation Results and 
Discussions 

We simulated sputtering process of hydrogen 
atoms and graphene sheets with the three 
incident-energy cases, namely Ei = 5 eV, Ei = 
15 eV and Ei = 100 eV. Each incident en- 
ergy corresponds to absorption process, reflec- 
tion process or penetration process of sputtering 
between a single graphene sheet and a hydrogen 
atom, respectively [2]. We show snapshots of 
graphene multilayer structure and injected hy- 
drogen atoms in the above three cases in Figs, 
n 01 and El respectively We also plot the radial 
distribution functions of carbon atoms for the 
first and the third sheets from the top in Figs. 
12 21 and |B1 Detail of each simulation result is 
shown in the foUowings subsection. 

3.1 The case of = 5 eV 

This situation corresponds to the absorption 
process between a hydrogen atom and a sin- 
gle graphene sheet [2]. It was intuitively ex- 
pected that incident hydrogen atoms destroy the 
graphene multilayer structure from the surface 
graphene sheet by breaking covalent bonds be- 
tween carbon atoms. The present simulation re- 
sults shows that our prediction is partially cor- 
rect and that the surface graphene, which is ex- 
posed to hydrogen atoms, is broken first. Con- 
trary to our prediction, however, the leading 
cause of destruction of graphene sheet is not that 
the incident hydrogen breaks the covalent bond 



between carbon atoms, but that the incident hy- 
drogen transfers its momentum to the carbon 
in the first graphenc sheet (Fig^a)). The first 
graphene sheet, which gained momentum from 
the incident hydrogen, undulates more deeply, as 
hydrogen atoms are injected. In course of time, 
amplitude of undulation becomes the same large 
as the distance of layers and the first graphene 
sheet contacts with the second one; then, the 
carbon in the first sheet interacts with the car- 
bon in the second sheet (Fig^^b)). Thus, the 
first and the second sheets deform and their 
honeycomb-structures are broken. The simula- 
tion in the next case of Ei = 15 eV makes it 
clear that breaking covalent bonds between car- 
bon by hydrogen do not play an important role 
in the above destruction mechanism. 

In Fig|2Ia), we plot the radial distribution 
functions g{r) among carbon atoms in the first 
graphene sheet for the three cases nH=l, nH=79 
and nH=159, where njj denotes the total number 
of injected hydrogen atoms. This figure shows 
that, as injected hydrogen atoms increases, the 
graphene structure of the first graphene sheet is 
breaking. On the other hand, Figl^b) shows 
that g{r) among carbon atoms in the third 
graphene sheet remains almost unchanged dur- 
ing hydrogen injection, which denotes that the 
third sheet is not destroyed by hydrogen. 

3.2 The case of Ei = 15 eV 

The second situation corresponds to the reflec- 
tion process of a hydrogen atom from a single 
graphene sheet [2]. Therefore, it is expected 
that it is difficult for hydrogen to connect with 
graphene sheet. This prediction is confirmed by 
Figl2a), which shows that the number of over- 
hanged hydrogen atoms is only four in spite of 
the fact that nn = 30. From this figure and 
Fig^a), deformation of the first graphene sheet 
is recognized in the case that nn — 30. By the 
above fact and the simulation result in the case 
that El = 5 eV, we reached the following conclu- 
sion: breaking the covalent bonds between car- 
bon atoms by hydrogen does not play an impor- 
tant role during destruction process of graphene 
structure, but momentum transfer from incident 
hydrogen to graphene causes to break graphene 
structure. 

Furthermore, we obtained new information 
about structures of hydrocarbon and carbon 
molecules which are yielded by sputtering. From 
FigEIb), it is found that almost all fragments 
of graphene sheets form chain-shaped molecules, 
and that yielded hydrocarbon molecules are 
composed of carbon chain and single hydrogen- 



atom. 

3.3 The case of Ej = 100 eV 

The last situation corresponds to the penetra- 
tion process of a hydrogen atom from a sin- 
gle graphene sheet [2]. Injected hydrogen pene- 
trates several graphene sheets with losing their 
kinetic energy, which is transferred to graphene 
sheets, until their energy becomes of the order 
of a few ten eV. After that, the hydrogen atoms 
begin to connect with carbon atoms (Fig[5fa)). 
From the standpoint of graphene sheets instead 
of hydrogen, only one side of graphene sheet is 
exposed to hydrogen injection in the cases that 
£"1=5 eV and i?i=15 eV. However, in the case 
of El = 100 eV, both sides of graphene sheets 
are attacked by hydrogen atoms. Therefore, 
graphene structure is destroyed rapidly in a wide 
range (FigEl^b) and (c)). Figure El shows that 
the third sheet is destroyed in the almost same 
timing as the first sheet. 

At the last stage of the simulation, all carbon 
molecules form chain-shaped structure. 

4 Conclusion 

Destruction process of graphene multilayer due 
to hydrogen injection is demonstrated by CMD 
with the modified REBO potential. From our 
simulation, we clarify the destruction process 
of graphene sheets as follows: (1) hydrogen 
atoms transfer their momentum to graphene 
sheet; (2) the graphene is undulated and touches 
the neighboring graphene sheets; (3) inter-sheet 
interaction of carbon atoms works and breaks 
graphene structure. Moreover, it is found that 
almost all fragments of graphene sheets form 
chain-shaped molecules, and that yielded hydro- 
carbon molecules are composed of carbon chain 
and single hydrogen-atom. 
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Fi gure 1: Snapshots of graphene multilayer struc- 
ture and injected hydrogen atoms in the case that Ei=5 
eV. (a) About fifteen hydrogen atoms connect with the 
graphene sheet in the case that the total number of the 
injected hydrogen atoms nu is 79. Each hydrogen atom 
and the surrounding carbon atoms compose the stable 
"overhanged" state; (b) njj has increased to 159 and the 
first graphene sheet contacts with the second one. The 
graphene structure is damaged. The increase in the num- 
ber of the overhanged hydrogen atoms is not as large as 
the increase in nh, i.e., 159 — 79 = 80. The third sheet 
remains almost unchanged. 
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Figure 2: Radial distribution functions g{r) among 
carbon atoms in the first (a) and the third (b) graphene 
sheets for the three cases nH=l, nH=79 and nH=159, 
where njj denotes the total number of injected hydrogen 
atoms. The incident hydrogen energy is 5 eV. 




Figure 3: Snapshots of graphene multilayer structure 
and injected hydrogen atoms in the case that i?i=15 
eV. (a) Four hydrogen atoms connect with the graphene 
sheet in the case that the total number of the injected hy- 
drogen atoms nil is 30. Each hydrogen atom and the sur- 
rounding carbon atoms compose the stable "overhanged" 
state; (b) njj has increased to 159 and the first graphene 
sheet contacts with the second one. Thus, the graphene 
structure is damaged. Fragments of graphene sheets and 
hydrogen atoms form chain-shaped molecules. The third 
sheet remains almost unchanged. 
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Figure 4: The radial distribution functions g{r) among 
carbon atoms in the first (a) and the third (b) graphene 
sheets for the three cases nH=l, nH=79 and nH=159. 
The incident hydrogen energy is 15 eV. 





Figure 5: Snapshots of graphene multilayer structure 
and injected hydrogen atoms in the case that i?j=100 
eV. (a) After injected hydrogen atoms penetrate sev- 
eral graphene sheets, hydrogen atoms contact with the 
graphene sheet in the case that the total number of the 
injected hydrogen atoms nji is 8. Graphene sheet is con- 
nected in both sides with hydrogen, (b) njj has increased 
to 30 and the upper five graphene sheet contact with 
other sheets; (c) graphene structure has already been 
broken to chain molecules when rth becomes 159. 
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Fi gure 6; Radial distribution functions g{r) among 
carbon atoms in the first (a) and the third (b) graphene 
sheets for the three cases nH=l, nH=79 and nH=159. 
The incident hydrogen energy is 100 eV. The function 
g{r) of the third sheet is similar to that of the first sheet. 



